Analysis of the atom-number correlation function in a few-atom trap 
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Stochastic properties of loading and loss mechanism in a few atom trap are analyzed. An ap- 
proximate formula is derived for the atom-number correlation function for the trapped atoms in the 
limit of reasonably small two-atom loss rate. Validity of the approximate formula is confirmed by 
numerical simulations. 
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I. INTRODUCTION 



Techniques for trapping small number of atoms in a 
microscopic volume have recently become an important 
tool for wide range of experiments in atomic physics and 
quantum optics such as cold collisions [3, atom metrol- 
ogy 2\ , cavityquantum electrodynamics Q and quantum 
information |^. 



There have been numerous measurements of trap load- 
ing and loss parameters of a magneto-optical trap (MOT) 
with many atoms, but they used indirect methods such 
as fitting to a model curve, dn/dt = B — Tn — J f3n^ 
|3, where n is the density of trapped atoms, T the loss 
rate due to collisions with background atoms, and (3 a co- 
efficient for two-atom collisional process among trapped 
atoms. Recently, several groups have trapped a small 
number of atoms in a MOT with a strong magnetic field 
gradient and could observe individual loading and loss 
events in real time. In this way, loading rate R, one- 
atom loss rate Li, _aiid two-atom loss rate L2 have been 
directly measured 



The atom-number correlation function has also been 
measured from the observed sequence of instantaneous 
atom number in a trap [8J. Since no analytic solution 
is known for the master equation in the limit of non- 
negligible two-atom loss, the study of the atom-number 
correlation function has been limited to the case of no 
two-atom loss, for which the correlation function does 
not provide any further information other than the one- 
atom loss rate. 

In the present work, we provide a comprehensive frame 
work for investigation of the atom-number correlation 
function based on the master equation (Chaps. II and 
III). An approximate formula has been derived for the 
correlation function in the limit of non-negligible two- 
atom loss (Chap. IV) and the validity of the approximate 
formula has been put to a test by numerical simulations 
(Chap. V). 
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II. DESCRIPTION OF MODEL 

The time dynamics of the atom number iV in a trap 
is governed by loading and loss processes. The loading 
process occurs at a certain rate R, called loading rate, 
which is determined by capture ability of the trap. The 
capture ability depends on such experimental conditions 
as laser intensity, laser beam size, laser-atom detuning 
and background source density of atoms but not on the 
number of atoms trapped already. We can assume R as 
a constant under fixed experimental condition. 

The loss rates, on the other hand, are affected by the 
number of atoms in the trap. One-atom loss occurs when 
one of the trapped atoms collides with a fast-moving 
background atom of different kind. One atom loss rate is 
thus linearly proportional to the number of atoms N al- 
ready present in the trap. We define Fi as the one-atom 
loss coefficient in such a way that the one atom loss rate 
Li is given by Li = FiA^. 

Two-atom loss process is due to the collision between 
two of the trapped atoms. There are several types of 
two-atom collision processes responsible for two-atom 
losses. They are ground-state hyperfine-changing colli- 
sions, ground-excited-state fine-structure changing colli- 
sion and ground-excited-state radiative escape. When- 
ever these collisions occur, the colliding two atoms can 
ain enough kinetic energy for escaping from the trap 
Therefore, the two-atom loss rate is proportional 
to the number of two-atom combinations for the atoms 
in the trap. In terms of T2 the two-atom loss coeffi- 
cient, the two-atom loss rate L2 is then given by L2 — 
T2N{N-l)/2. 

The loading and loss processes occur randomly. Due to 
random nature of these processes, it is more convenient to 
treat the problem in terms of the atom number distribu- 
tion function Pjv than to deal with the time variation of 
the instantaneous atom number N itself. In order to un- 
derstand the connection between time evolution of P^r {t) 
and that of N , let us suppose that we turn on the load- 
ing process at time t — Q and observe the atom number, 
initially zero, afterwards. The atom number will change 
randomly as depicted in Fig. Q] but the center of fiuctu- 
ation will increase to a steady-state value determined by 
the balance between the loading and loss processes. 

If we repeat the observation infinitely many times, each 
observation will give a different time sequence in detail. 
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FIG. 1: (a) Ensemble picture for the loading process turned 
on at f = 0. For each member of the ensemble, the atom num- 
ber constantly fluctuates due to random loading and loss pro- 
cesses, (b) Ensemble averaged atom number in time, (c) En- 
semble averaged atom number distribution at various times. 



However, if we average all of the observed time sequences 
with the initial starting time ahgned, we obtain a distri- 
bution of atom numbers at any time t, corresponding 
to PAr(i), and a sequence of averaged atom number as 
a function of time. Alternatively, we can replace the in- 
finitely many observations with an ensemble of identically 
prepared traps. The sequence of averaged atom number 
is then an ensemble average of the atom number (N) . 

It is shown that the time evolution of the ensemble- 
averaged distribution Pjy {t) can be described the follow- 
ing master equation: 



A^+2 



A^+l 



N 



P 
P 



N-\ 



R ' 


r,(A/+2)j 


1 

1 r2(A^+2)(A^+l)/2 


R 


r,(A^+i)i 


1 1 

1 ^T,(N+\\Nn 


R 


''"i. 


1 1 

, .YMN-m 

▼ 




R^ Y,(N-\)\ 


1 



N-\ 



FIG. 2: Connections among probabilities Pjv's. Loading pro- 
cess TV — 1 — > A'', one-atom loss process A^ + 1 — + A'^ and 
two-atom loss process N + 2 ^ N result in an increase in Pjv 
whereas loading process A'' -^ A'^ -I- 1, one-atom loss process 
N —> N — 1 and two-atom loss process N —> N — 2 cause a 
reduction in Pjy. 



III. WITHOUT TWO-ATOM LOSS TERMS 

If a density of atoms in a trap is low enough, the col- 
lisions among the trapped atoms can be neglected. If we 
assume r2 = 0, Eq. (Q) is simplified as 



dP, 



N 



dt 



= RPn-i+Ti{N + 1)Pn+i - {R + TiN)Pn . (2) 



This equation corresponds to the well-known birth-death 
model 9]. In the steady state, we have dP^/dt ~ for 
all N and the solution is a Poisson distribution given by 



^- = M^ 



-N^N 



(3) 



for all integers A^ > with N = R/Ti, the mean atom 
number. 

A rate equation for the ensemble-averaged atom num- 
ber (A^) in the trap can be derived from the Eq. JSJI. 
Multiplication of A^ to the both sides of Eq. (0) followed 
by summation over N gives a differential rate equation 
for (A^): 



dP. 



N 



dt 



RPn^i+Ti{N+1)Pn+i 
+ ^T2iN + 2){N + l)PN+2 
-{R + r,N + ^N{N - 1)T2}Pn 



(1) 



for A^ = 0, 1, 2, . . . with a convention of P_i = 0. Con- 
nections among probabilities Pat's are depicted in Fig. 
11 



dt 



(N)=R-Ti{N) 



(4) 



where the ensemble average is formally defined as 

oo 

(/(TV))^^Pa./(TV). (5) 
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The formal solution of Eq. ^ in terms of an initial atom 
number A^o is given by 



(Nit)) = -J 
i 1 
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(6) 



In the steady state, d{N)/dt = and thus we get (N) = 
N = R/Ti from Eq. Q. Here the upper bar indicates 
a time average in the steady state. From the properties 
of the Poisson distribution we get the following relation 
between the variance and the mean. 
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(7) 



In the steady state, the correlation function of the atom 
number is defined as follows. 



CiT)^{N{t)N{t + T))t 



(8) 



where r is a time delay and the notation (. . .)t represents 
a time average. Although the ensemble-averaged atom 
number {N) does not change in the steady state, equal 
to the mean atom number N, the atom number A'^ itself 
is continuously fluctuating around its mean. 

We can replace the time average above with an ensem- 
ble average. Let us denote the atom number at time t 
in the ith ensemble member as N^'^'{t). Then we can 
rewrite the correlation function as 



C(r) = ^iV(')(0)A«(r) 



(9) 



In the second line in Eq. H12() . we used Eq. Q, which is 
independent of the initial condition. By using the con- 
ditions C(0) = (N^)t =W and N = R/Ti, we obtain 



dC{T) 



dr 
and 



T = 



RN - TiN^ ^RN ~ Ti{N^ + N) 



C{t) =N^ + (t'^ 



and therefore, the normalized correlation function is 
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-R 

(13) 

(14) 

3 

(15) 



The correlation function shows a decay behavior from 
1 + 1/iV to 1 with a characteristic correlation time tc = 
l/Fi. It is interesting to note that the atom-number cor- 
relation function exhibits bunching, i.e., C(0) > C{t) for 
T > 0, when the atom number distribution is Poissonian 
whereas a Poisson distribution for photon number does 
not necessarily means bunching for a photon number cor- 
relation function. 



where the summation over i is performed over all mem- 
bers of the ensemble. We can regroup the ensemble into 
sub-ensembles in such a way that the members in each 
sub-ensemble have a common initial atom number Nq. 
Since A''^*) values are distributed according to P^ of Eq. 
(3), we can rewrite the above equation as 



C(r)-^AoPAr„ 

JVo 



J2n('\. 



N('"l(0) = No 



(10) 



where the summation over i' represents a summation over 
each sub-ensemble. All sub-ensembles are statistically 
identical, described by the same master equation, Eq. (2). 
Now, we can see that the quantity in [. . .] is nothing but 
the ensemble-averaged atom number (A^) at time t = t 
when its initial value is A'o at time t = 0. The correlation 
function is then simplified as 



C{t)=Y,N^P^AN(t))\ 



{N{0})=No 



(11) 



No 



A differential equation for the correlation function can 
be obtained by taking a derivative of Eq. Hll|) with re- 
spect to T. 



dr ^'--^"Ut 

No 

oo 

J2NoPn„ {i?-ri (A^(r))|^^(o„^^J 



No 

= RN-TiCir) 



(12) 



IV. WITH TWO-ATOM LOSS TERMS 

If the density of atoms in a trap becomes higher, the 
two-atom loss terms are no longer negligible, and thus 
the full version of the master equation, Eq. (Q, should 
be considered. However, the full master equation cannot 
be solved analytically. Instead, we rely on numerical so- 
lutions, either by Monte-Carlo simulation or an iteration 
method for given parameter values. Such numerical stud- 
ies are discussed elsewhere [i3- In this work, we focus 
on approximate solutions to the master equation. 

The rate equation for (A^) with the two-atom loss terms 
can be derived from the master equation as 



dt 



(N) =R- Ti{N) - r2{N{N - 1)) 



(16) 



In the steady state, we let d{N)/dt = 0. Differently 
from the case of excluding two-atom loss terms, however, 
we do not get any information about A^ other than the 
following relation between the variance and the mean. 
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(17) 



Following the same line of reasoning from Eq. © to 
Eq. H12I) with the rate equation Eq. 1)16(1 . we can obtain 
a differential equation for the correlation function in the 
presence of non-negligible two-atom loss terms as follows. 



dCJT) 
dr 



RN+{r2-Ti)C{T)~T2{N{t)N\t + T)) (18) 



When r = 
comes 

dC{T) 



dr 



0, the slope of the correlation function be- 



= -RN + {r2 -Ti)N^ -T^N^ , (19) 



where N^ — {N^)t- The relation of N^ with the lower 
moments is obtained from the rate equation for (iV^) in 
the steady state. 

(N^) = = R{2N+l)-Ti(2W-N)-2r2{W-2W+N) 

(20) 
After some lengthy algebra we finally obtain 



dt 



dC{T) 



dr 



2 2 



(21) 



If r2 = 0, Eq. (EU reduces to Eq. (P|l . 

The differential equation Eq. (|18|l cannot be solved ex- 
actly because of the complexity in the last term. How- 
ever, in the limit of reasonably small two-atom loss, 
r2 <C R,Ti, the atom number distribution is approxi- 
mately a Poisson distribution with a variance a^ « N. 
This is the most significant approximation in our analy- 
sis. The validity of this approximation will be discussed 
in the next section. In this limit, the mean value N can 
be obtained from Eq. (|17|l . 
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and thus 



N : 



-ri + Vrf + 4lr2 



2R 



2T2 



N' 



(23) 



4i?r2 



Ti + r 

Equation H23|l reduces to Eq. 



where F' = yTf 
(7) when Ta -^ 0. 

Our second approximation is that the correlation func- 
tion is in the same form as that of Eq. (|14|l except for 
the decay constant. This is a reasonable approximation 
since the correlation function always starts from A^^ and 
goes to N'^ at infinity. Therefore, we assume 
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(24) 



where Feg is an effective decay constant to be deter- 
mined. It can be determined by Eq. H21I) . 
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and thus 



'='' - 2iV 2 



El + sr 



(26) 



The normalized correlation function is then given by 
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Using the approximation a^ 

1 



C(r) « 1 + 
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iV, we obtain 
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FIG. 3: Comparison of atom-number correlation functions by 
numerical simulation (solid curve) and by approximate for- 
mula (dotted curve). Parameters are _R=0.08, ri=0.04 and 
the two-atom collision constant r2 is (a) 0.004, (b) 0.04 and 
(c) 0.4, respectively. A dashed curve represents a correlation 
function given by Eq. 1151 neglecting two-atom loss. Grey 
curves are the fit given by the approximate formula with a^ 
treated as a fitting parameter. The resulting a^ values are 
0.94iV, 0.78iV and 0.50iV for (a), (b) and (c), respectively. 



COMPARISON WITH SIMULATION 
RESULTS 



In order to check the validity of our approximate for- 
mulas Eqs. (|23|l and (|28|) . we solve the master equa- 
tion numerically and compare the results with the ap- 
proximate formulas in Fig. |31 Under the condition of 
F2 ^ .R, Fi, Eq. (|28|l agrees well with the numerical re- 
sults as shown in Fig.|3Ja). Our approximation is valid in 
this limit. As F2 is increased beyond this limit, however, 
the approximation starts to deviate from the simulation 



result as seen in Fig. [Sfb) , where we have r2 ~ -R, Fi . 
The approximation is still acceptable even under this 
condition although not as good as in Fig. Ofa). The 
deviation is severe in Fig. ^c), under the condition of 

r2>i?,ri. 

The failure of the approximation is mostly due to the 
fact that the atom number distribution is assumed to be 
a Poissonian. The atom number distribution function is 
not well approximated by a Poisson distribution when r2 
becomes comparable to and larger than R and Fi , so the 
assumption, N ^ a"^, fails. In general, the variance a^ 
is smaller than the mean value N. Consequently, when 
T — 0, the approximate result is always larger than the 
true value. 
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CiO) (29) 



If we treat a^ in Eq. (|27|l as a fitting parameter, we obtain 



better agreement between the numerical results and the 
approximate formula even for r2 ^ Fi, R. For example, 
the grey curves in Fig. 01 are the fit given by Eq. l(?7|l . 
agreeing well with the numerical results. 



VI. CONCLUSIONS 

We have derived approximate formulas for the mean 
atom number and the atom-number correlation function 
in the limit of reasonably small two-atom collision com- 
pared to the one-atom collision and the loading rates. 
The validity of the approximate formulas was confirmed 
by comparing them with numerical solutions. 
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